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Narrow-band detection of the Raman water vapor spectrum using the lidar technique introduces a 
concern over the temperature dependence of the Raman spectrum. Various groups have addressed this 
issue either by trying to minimize the temperature dependence to the point where it can be ignored or by 
correcting for whatever degree of temperature dependence exists. The traditional technique for perform- 
ing either of these entails accurately measuring both the laser output wavelength and the water vapor 
spectral passband with combined uncertainty of approximately 0.01 nm. However, uncertainty in inter- 
ference filter center wavelengths and laser output wavelengths can be this large or larger. These com- 
bined uncertainties translate into uncertainties in the magnitude of the temperature dependence of the 
Raman lidar water vapor measurement of 3% or more. We present here an alternate approach for 
accurately determining the temperature dependence of the Raman lidar water vapor measurement. This 
alternate approach entails acquiring sequential atmospheric profiles using the lidar while scanning the 
channel passband across portions of the Raman water vapor Q-branch. This scanning is accomplished 
either by tilt-tuning an interference filter or by scanning the output of a spectrometer. Through this 
process a peak in the transmitted intensity can be discerned in a manner that defines the spectral 
location of the channel passband with respect to the laser output wavelength to much higher accuracy 
than that achieved with standard laboratory techniques. Given the peak of the water vapor signal 
intensity curve, determined using the techniques described here, and an approximate knowledge of 
atmospheric temperature, the temperature dependence of a given Raman lidar profile can be determined 
with accuracy of 0.5% or better. A Mathematica notebook that demonstrates the calculations used here is 
available from the lead author. © 2013 Optical Society of America 
OCIS codes: (010.3640) Lidar; (300.6450) Spectroscopy, Raman. 

http://dx.doi.org/10.1364/AO.52.005376 


1. Introduction 

Raman lidar is a well-established technique for 
probing water vapor in the atmosphere dating from 
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the first nighttime measurements of Cooney [1] and 
Melfi [2] to the current capability of day and night 
measurements offered by numerous lidar sytems 
[3-8]. Although very useful Raman lidar measure- 
ments of water vapor and other species have been 
made with the solar blind technique [9] , the narrow- 
band, narrow field of view (NBNFOV) approach, 
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which makes use of high UV laser energy, narrow 
spectral and spatial filtering and a combination of 
analog and photon counting data acquisition [3,4], 
has become the most popular Raman lidar technique 
for measuring water vapor. The reason for the popu- 
larity of the NBNFOV technique is that it permits 
useful daytime measurements to be performed with- 
out limiting the upper range of nighttime measure- 
ments as is the case when using the solar blind 
technique. 

It was first noted by Sherlock et al. [10], however, 
that the use of narrow band detection of the water 
vapor spectrum introduces a concern about the 
temperature dependence of the Raman lidar water 
vapor measurement. That work was based on a 
Raman scattering model that did not include aniso- 
tropic scattering. Fundamental theoretical work in 
the spectroscopy of the Raman OH stretch region, 
performed by Avila et al. [ 11 ] and that built upon 
the work of Murphy [12,13] and Fernandez-Sanchez 
[14,15], included anisotropic scattering explicitly. 
This advance in Raman scattering modeling permit- 
ted Whiteman [16,17] to demonstrate an accurate 
correction technique for this temperature depend- 
ence of water vapor using a modified version of the 
lidar equation that explicitly contained the temper- 
ature dependence as a multiplier of the traditional 
lidar equation. It became clear from that work, 
however, that laser and interference filter character- 
istics had to be known with quite high accuracy to 
optimally implement the correction. Whiteman et al. 
applied this formulation to the correction of a 
wet bias observed during the International H 2 0 
Project (IHOP) field campaign [4] and Atmospheric 
Infrared Sounder (AIRS) validation measurements 
[ 18 ] that significantly improved the upper tropo- 
spheric measurements of their scanning Raman 
lidar. Application of the technique, however, was per- 
formed using high resolution wavemeter measure- 
ments of the output spectrum of the laser, as 
well as mercury lamp calibrated Fourier transform 
spectrometer (FTS) and spectrophotometer measure- 
ments of the interference filter. It was also necessary 
to assess any deviations from pure normal position- 
ing of the interference filter. 

Given that discrepancies between manufacturer’s 
filter transmission curves and high resolution FTS 
measurements have shown differences typically of 

0.01-0.02 nm, but at times as large as 0.05 nm [18], 
implementation of an accurate solution for this tem- 
perature dependence requires very careful labora- 
tory work whether the desire is either to minimize 
the temperature dependence or to correct for it. 

Because of the difficulty of determining the tem- 
perature dependence of Raman lidar measurements 
of water vapor using these previously described 
techniques, a new approach was conceived and 
implemented in 2007 as part of the NASA/GSFC 
ALVICE mobile lidar system [5]. That approach 
entails tilt-tuning [ 19 ] interference filters for both 
the Raman water vapor and nitrogen channels while 


the filters remain in the optical train and acquiring 
a sequential set of atmospheric profiles. Although 
the discussion here will focus mainly on the use of 
interference filters, the same approach can be imple- 
mented using a spectrometer-based system by scan- 
ning the output of the spectrometer. Also, in this 
article we consider only the temperature dependence 
of the Raman water vapor measurement. The meas- 
urement of the water vapor mixing ratio from Raman 
lidar also requires measuring the Raman nitrogen 
signal, which possesses some temperature depend- 
ence as well. But, provided that narrow-band filters 
with widths in the range of 0.2-0. 3 nm are used, the 
temperature dependence of the Raman nitrogen 
measurement is less than 0.5% over the range of 
normal atmospheric temperatures [ 16 ] and thus 
will not be considered here. Note, however, that the 
use of wide-band filters in the measurement of the 
Raman nitrogen signal can increase the temperature 
dependence of the measurement into the 1%— 2% 
range [16]. We will describe the new approach pro- 
posed here in more detail after reviewing the 
temperature-dependent lidar equation [16,17], 

2. Temperature-Dependent Lidar Equation 

The idea behind this equation, which was motivated 
by the work of Sherlock et al. [10], is to include 
the temperature dependence of the lidar equation 
into a separate, multiplicative factor that can be as- 
sessed separately from the elements of the normal, 
temperature-independent, lidar equation. That fac- 
tor is referred to as F H below. 

The temperature-dependent lidar equation can be 
written 

P(r . X L , X H , AA H ) = Po 2 i](/i H )0(r )F H ( T)/i^ (X L , r) 

2 r z 

x exp J [a(A L . r') + a(A H . r'jjdr' j 
+ B(r,Ajj, A/hf), (1) 

where 

1. P(r, A L ,A H , AA h ) is the return signal power 
measured at range r due to laser emission at A L using 
a bandpass centered at A H with bandpass width A A H . 

2. P 0 is the laser power. 

3. ct/2 is one half of the length of the volume 
illuminated by the laser beam with r being the laser 
pulse temporal width and c the speed of light. 

4. A/r 2 is the solid angle subtended by the 
receiver as viewed from the scattering volume. 

5. r](A H ) is the optical efficiency of the water 
vapor channel evaluated at the bandpass center 
wavelength (CWL) of Ah- 

6. F h (T) is the factor in the lidar equation that 
carries all the temperature dependence of the Raman 
water vapor measurement and is defined as 
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f^ H nW° H {X,n,T)/dQ.&{X) 

rj{X H )da H {n) / dQ. 

where do H (n)/d£l is the full Raman backscatter cross 
section, here restricted to the Q-branch due to the 
width of the filter considered, when excited at the la- 
ser wavelength. da H (X. n. T)/d£i is the temperature- 
dependent Raman backscatter cross section where 
we note that the individual lines within the Q-branch 
are temperature dependent but the integral over the 
full Q-branch, given by do H (n) / dQ is essentially 
temperature independent. >/(/) is the optical effi- 
ciency of the channel passband as a function of wave- 
length, X, within the intervalAl# that is defined by 
the filter passband. In the numerical work to follow, 
the integral in the numerator is implemented as a 
simple summation [ 20 ] over the lines tabulated in 
Avila et al. [ 21 ] F H provides the fraction of the full 
Raman Q-branch cross section that is transmitted 
by the Raman lidar channel passband (e.g., interfer- 
ence filter) given the atmospheric temperature of the 
scattering volume and the transmission characteris- 
tics of the passband. 

7. fi^(A L .r) — N H {r)do H {n)/dil is the Raman 
backscatter cross section due to excitation, at the 
laser wavelength, of N H concentration of water vapor 
molecules at range r. 

8. exp (-/q[«(/Il, r ') + a(A H , r')dr']) is the round- 
trip transmission to the scattering volume and back 
and involves the integral of the extinction, a, at the 
laser wavelength on the outbound path and at the 
Raman shifted wavelength on the return path. 

9. B(r, A h ,AA h ) is the background power due to 
skylight and detector noise. 

It is the F h term defined in number 6 above that 
will receive much of the focus of the following sections. 

3. Simulations of the Raman Water Vapor Spectrum 

The temperature-dependent correction discussed 
here has been applied in different contexts. It has 
been used to improve the accuracy of water vapor 
data acquired during field campaigns, such as the 
International H 2 0 Project [4], the AIRS Water Vapor 
Experiment-Ground [18], the Water Vapor Experi- 
ment Satellite/Sondes experiment [6], and for ex- 
tended meteorological observations in Japan [ 22 ] . It 
has also been used as part of the implementation of 
an independent calibration of a Raman water vapor 
lidar [23,24]. All of these applications have been 
based on the tripled Nd:YAG laser (~354.7 nm) and 
that is also the focus of this article. 

It was noted before [4] that, for a laser output 
wavelength of 354.7 nm, tuning a 0.25 nm wide inter- 
ference filter to a wavelength ~407.45 nm (Raman 
shift of ~3649 cm -1 ) would minimize the temperature 
dependence of the Raman signal from water vapor 
although the exact temperature-independent location 
in the spectrum is a function of the filter shape. Also, 
with typical uncertainties in the CWL of the interfer- 
ence filter in the range of 0.01-0.05 nm, it becomes 


difficult to time an interference filter with confidence 
to either minimize the temperature dependence or to 
maximize the signal intensity, where we note that sig- 
nal intensity is sacrificed if the filter is tuned to the 
nearly temperature insensitive point of ~3649 cm -1 . 
We will illustrate this after first demonstrating 
calculations of the Raman Q-branch using the data 
provided by Avila et al. [20]. 

Avila et al. 2004 [ 20 ] is an update to the 1999 work 
of the same Spanish group. The main advances in the 
2004 work are due to improved ro-vibrational wave- 
functions and are seen at temperatures exceeding 
400 K and regions outside of the Q-branch. To illus- 
trate the use of the published transitions, the Raman 
Q-branch was simulated at 297 Kwith 0.14 cm -1 res- 
olution and is shown in Fig. 1. Two filters with 
approximate widths of 0.25 nm (the first from Fig. 3 
of Whiteman et al. [4] and here tuned to 3652 cm -1 ” 
the second to be described later) are also shown for 
comparison purposes. The first of these will be the 
spectral shape used for most of the calculations to 
follow. The Q-branch spectrum shown in Fig. 1 
compares favorably with Fig. 4(b) of Avila et al. 
2004 [ 20 ] and is essentially indistinguishable from 
the spectrum calculated using the results from Avila 
et al. 1999 [ 11 ] (Fig. 1 in [16]). The calculation of the 
full Raman Q-branch cross section performed here 
also agrees to within 0.001% of the values shown 
in Table 3 of Avila et al. 2004 [ 20 ] when using the 
same emission wavenumber and partition function 
values [25]. It should be noted that, based on theo- 
retical considerations and confirmed through calcu- 
lations using these 2004 data, the full Q-branch 
intensity is essentially independent of temperature 
over the range of temperatures of interest in Raman 
lidar measurements of the normal atmosphere. 

We now focus on calculating the F h term shown in 
Eq. (1) as this term carries all the temperature 
dependence of the Raman lidar equation. As 
mentioned, F H can be considered as the effective 
Raman Q-branch cross section that results from 


Q-branch Filterl Filter2 



Fig. 1. Calculation of the Raman water vapor Q-branch simu- 
lated with 0.14 cm -1 resolution at 297 K in order to replicate 
Fig. 4(b) in Avila et al. 2004 [20] . Two filter shapes are shown 
corresponding to the filters used in experiments performed in 
2007 and 2011, respectively, and are described later in the text. 
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Filter Center (cm' 1 ) 


Fig. 2. Calculation of F H for Filterl shown in Fig. 1. but with 
the center at different positions in the water vapor Q-branch. 
Calculations are performed over the range of 200-300 K. The peak 
in the F H curve shifts approximately 1 cm -1 over the 100 K range 
considered. Note that a filter located at a shift of 3649.5 cm -1 lim- 
its the temperature sensitivity to approximately 1% over the range 
of temperatures considered, but that the transmitted intensity for 
a filter positioned at this location is ~4-12% less than if the filter 
were positioned at the peak of the F H curve. 

the convolution of the individual temperature- 
dependent lines in the Q-branch spectrum and the 
spectral passband in use. Figure 2 shows calcula- 
tions of F h using the shape of Filterl shown in 
Fig. 1 where that shape is centered over a range of 
3648-3656 cm -1 , simulating a tilt-tuning experi- 
ment. Temperatures between 200 and 300 K are 
considered given that this range covers most temper- 
atures encountered in the troposphere and lower 
stratosphere. The largest effective cross sections are 
achieved at CWL positions between approximately 
3652 and 3653 cm -1 depending on temperature. This 
amounts to a total shift of approximately 1 cm -1 in 
the peak of the F H curve over the 100 K temperature 
range considered. The figure indicates that at an 
atmospheric temperature of 300 K, using Filterl 


• Calc_ 285 ■ 5/21/2007 
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Fig. 3. (Blue) Theoretical calculation of F H (T) using Filterl 
shown in Fig. _1 for a temperature of 285 K (Media 1). (Magenta) 
Tilt-tuning lidar measurements from the night of 21 May 2007 
of transmitted intensity in the lidar water vapor channel as a 
function of tilt angle normalized to the theoretical F h (T) curve. 


• C'alc_ 278 K ■ 10/31/2011 a ♦ 10/31/20111) 



Fig. 4. (Blue) Theoretical calculation of Fh (T) using Filter2 shown 
in Fig. 1_ for a temperature of 278 K. (Magenta and Brown) Lidar- 
derived results from two similar tuning experiments performed on 
31 October 2011 using the tilt module described in the text. 

shown in Fig. 1 but considering it to have 100% 
peak transmission, almost 86% of the total Raman 
Q-branch cross section would be transmitted. Con- 
sidering a 200 K temperature, the same filter posi- 
tioned at 3653 cm -1 would transmit approximately 
92% of the Raman Q-branch, an increase of more 
approximately 7%. The plot also indicates that, for 
this filter, the transmitted intensity of the filter 
positioned at 3649.5 cm -1 would change by approxi- 
mately 1% over the 100 K range considered. Position- 
ing the filter passband at this location is the 
approach that has been taken by some groups [8] . 
But there are several reasons to prefer positioning 
the passband to the peak of the F H curve. 

As shown in Fig. 2, depending on temperature, 
positioning the filter to the nearly temperature- 
insensitive point of 3649.5 cm -1 sacrifices between 
4%— 12% of the transmitted intensity versus position- 
ing the filter at the peak of the F H curve. Positioning 
the filter at the peak of F H introduces additional 
temperature dependence when compared to position- 
ing the filter at ~3649.5 cm -1 but, as will be shown, 
performing an accurate correction for this tempera- 
ture dependence requires only an approximate 
knowledge of the atmospheric temperature profile. 
Finally, positioning the filter to 3649.5 cm -1 , which 
is along the shoulder of the F H curve, is more difficult 
experimentally than positioning the filter to the 
peak of the F H curve. 

4. Tilt-Tuning to the Maximum in the F h Curve 

The procedure suggested here is to acquire succes- 
sive lidar profiles through sequential tilt -tuning of 
the filter (or scanning the output of a spectrometer) 
and comparing these results with calculated Fh 
curves as shown in Fig. 2. This is the procedure that 
we have followed since first assessing the tempera- 
ture correction for the IHOP [4] and AWEX-G [ 18 ] 
field campaigns, held in 2002 and 2003, respectively, 
that were mentioned earlier. During the early imple- 
mentation of this technique, we used wedges of 
varying thicknesses to implement different angular 
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tilts of a mounted interference filter. At each angle, a 
separate lidar profile was acquired. This procedure 
required dismounting the filter between each lidar 
profile, which was a time consuming process. None- 
theless, under stable atmospheric conditions this 
approach provided reasonable results. An example 
of this is shown in Fig. 3. The theoretical F H curve, 
calculated using Filterl shown in Fig. 1 for a temper- 
ature of 285 K (see Fig. 3 multimedia file Media 1 for 
a pdf of the Mathematica notebook used to calculate 
this), is shown along with the results of a tilt-tuning 
lidar experiment performed on the night of 21 May 
2007. The temperature of 285 K was chosen to reflect 
the approximate temperature of the lower tropo- 
sphere during the time of the lidar measurements 
on this evening. To perform the lidar measurements, 
the water vapor filter was sequentially tilted in 1 
degree increments by using wedges of appropriate 
thicknesses. At each angle, water vapor profiles were 
acquired under stable nighttime conditions. To cre- 
ate the lidar derived values shown in Fig. 3, the sum 
of the counts in the water vapor channel between the 
altitude ranges of approximately 0.5— 2.0 km was 
formed for each profile. A best fit of these intensity 
values to the F H values shown was done to find the 
correspondence between tilt angle and Raman shift 
as shown in the figure. As the figure indicates, a 1 
degree tilt in the filter corresponded to a change in 
filter spectral location of approximately 1.5 cm -1 . 
The agreement between the calculated curve and the 
2007 measurements is quite good particularly in the 
shoulders. The close agreement of the shoulders 
permits the center of the filter to be determined with 
an estimated uncertainty of less than 0.2 cm -1 . 

Since 2010, we have used specially designed filter 
modules for both the water vapor and nitrogen chan- 
nels in the ALVICE Raman lidar. These modules 
permit the tilt-tuning experiment to be quickly and 
repeatedly implemented. The filter modules also per- 
mit accurate changes in tilt angle of less than 0.1 
degrees improving on the ability to assess the center 
location of the filter over the 2007 experiment shown 
in Fig. 3. Two experiments performed on the night of 
31 October 2011 using the tilt module are shown in 
Fig. 4 where the F H curve shown was calculated at 
the temperature of 278 K, the mean temperature 
on this night of the atmospheric layer used for the 
lidar measurements. Note that Filter2 from Fig. 1 
was used for these 2011 experiments. Most of the 
differences in the F H curves between Figs. 3 and 4 
are due to the difference in filter shape as the 
atmospheric temperature difference of 7 K does not 
produce a large change in the shape of the F H curve 
as can be discerned from Fig. 2. Performing the tilt- 
tuning experiment in a short period of time in gen- 
eral decreases the influence of changing atmospheric 
water vapor amounts on the results. The ability to 
acquire data in a shorter period of time, along with 
the higher angular resolution of the 2011 measure- 
ments, helps to explain the better correspondence 
between the theoretical curve and the measurements 


taken in 2011 versus those in 2007 shown in Fig. 3. 
This better agreement of curves reduces the uncer- 
tainty in the assignment of the center position of 
the filter to an estimated 0.1 cm -1 . 

Knowledge of the atmospheric temperature is 
needed to determine the spectral location of the in- 
terference filter through the experiments illustrated 
in Figs. 3 and 4. Note, however, that the peaks of 
the intensity curves shown in Fig. 2 change by only 
~1 cm -1 over the 100 K range considered. Therefore, 
only an approximate knowledge of the temperature 
of the atmospheric layer used for the summation of 
the lidar measurements is needed in order to deter- 
mine, with high accuracy, the appropriate F H curve 
to use for fitting the lidar-derived data. This is 
summarized in Table 1 where a 25 K change in the 
temperature of the atmospheric layer used for lidar 
measurements of signal intensity implies approxi- 
mately a 0.7 cm -1 change in the location of the 
peak in the F H curve. Under most experimental con- 
ditions the knowledge of the mean temperature of 
the atmospheric layer probed by the lidar can be de- 
termined to much better than 25 K. For example, if 
the mean temperature of the atmospheric layer 
probed is known to within 3 K, the uncertainty in 
the theoretical position of the peak of the F H curve 
in comparisons, such as shown in Figs 3 and 4 can be 
limited to less than 0.1 cm -1 which corresponds to 
approximately 0.001 nm when considering Raman 
lidar measurements of water vapor using the tripled 
Nd:YAG laser (~354.7 nm). 

5. Application of the Temperature Correction 

The temperature correction is determined by evalu- 
ating F H as a function of altitude and including F H as 
a multiplier in the lidar equation as shown in Eq. (1). 
Given the ability to evaluate F H for a single temper- 
ature, it is straightforward to evaluate F H for the 
entire lidar profile. But given that F h changes by 
6%— 9% over a 100 K range for the filter configura- 
tions considered in Table 1, it can be seen that the 
F h correction is not strongly a function of tempera- 
ture. To illustrate this point, consider Fig. 5. 

Figure 5 presents F H versus altitude for the situa- 
tion of Filterl shown in Fig. 1 centered at 3652 cm -1 , 
a situation similar to that achieved in the tilt-tuning 
calibration of the Raman lidar in the ALVICE system 


Table 1 Considering the 0.25 nm Filter Shown in Fig. 1 , the Peak in F„ as 
a Function of Filter Spectral Location; Percentage Change in Fh over the 
100 K Range of 200-300 K is Given 


Peak of F h Curve as a F unction of Atmospheric Layer and Percent 
Change in Fjj Over Range of Atmospheric Temperatures 

Atmospheric 

F h Peak Intensity 

Change (%) in F H 

Temperature (K) 

(cm -1 ) 

from 200 to 300 K 

200 

3653.2 

8.9 

225 

3652.9 

8.4 

250 

3652.6 

7.2 

275 

3652.3 

6.7 

300 

3652 

6.1 
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Fig. 5. F h versus altitude for the U.S. Standard Atmosphere [ 26 ] 
plus and minus 5 K using Filterl shown in Fig. 1. Note that to 
constrain the uncertainty in the temperature correction to less 
than 0.4%, the atmospheric temperature at a given altitude need 
not be known to better than approximately ±5 K 

when the focus is on lower tropospheric measure- 
ments. The U. S. Standard Atmosphere [ 26 ] is used 
for the temperature profile and deviations of ±5 K 
are considered at each altitude. Over the range of 
temperatures in the U. S. Standard Atmosphere from 
the surface to the tropopause, F H changes from 
approximately 0.863 to 0.902 for a total change of 
4.3%. For temperature changes of ±5 K considered 
at each altitude, the values of F H deviate by less than 
0.2%— 0.4%. Therefore, the figure illustrates that, 
provided the temperature is known to within ±5 K 
as a function of altitude, the uncertainty in the F h 
value due to uncertainties in atmospheric tempera- 
ture can be constrained to less than 0.4%. 

6. Discussion 

Performing the lidar tilt-tuning experiment and per- 
forming a best-fit of the lidar-derived values to the 
theoretically determined ones as shown in Figs. 3 
and 4 must be done with care. Note that the entire 
range of F H intensity values plotted in the figures 
is approximately 10%. Therefore, any changes in 
atmospheric water vapor amounts or laser output 
power during the time of the tilt-tuning experiment 
should be determined to much better than 10%. A 
sequence of lidar profile measurements of water 
vapor mixing ratio in the nominal configuration, 
taken both before and after the tilt-tuning experi- 
ment, can be used to assess the stability of the water 
vapor field during the experiment. Correlative mea- 
surements of total column water can also help to 
assess changes in water vapor amount during the ex- 
periment. Performing the experiment quickly, as is 
aided by the ability to tilt-tune the filter in place, also 
helps to minimize changes in water vapor amount. 
Assuming insignificant changes in the atmospheric 
density and extinction profiles during the experi- 
ment, changes in laser power are compensated for in 
the lidar measurements used to determine the mix- 
ing ratio. Note, however, that for the measurements 


shown in Figs. 3 and 4 only the water vapor signal 
intensity was monitored and not the mixing ratio. 
The stability of the measurements was such, how- 
ever, that no adjustments were needed for changes 
in water vapor concentration or in laser power. 

As mentioned earlier, our experience at GSFC is 
that wavelength uncertainties in manufacturer’s 
transmission curves are typically 0.01-0.02 nm. Also, 
the laser output wavelength changes as the atmos- 
pheric index of refraction changes due to changes 
in temperature and pressure [27]. One of the conse- 
quences of this is that there is a small change in 
wavelength as a function of pressure at a given eleva- 
tion even if the temperature is held constant. If the 
elevation at which the laser is operated changes, 
larger changes in wavelength can occur. For example, 
the tripled Nd:YAG laser operated at an elevation 
of 2.5 km will have an output wavelength approxi- 
mately 0.027 nm longer than if operated at sea 
level if only pressure changes are considered. Such 
a change in wavelength, if uncompensated for, and 
assuming a 0.25 nm filter can alter the water vapor 
mixing ratio calibration value of a Raman lidar by 
2%— 4% depending on the CWL of the filter. Because 
of these changes in laser output wavelength as a func- 
tion of atmospheric refractive index, careful estimates 
of the actual laser output wavelength include uncer- 
tainties in the range of 0.01-0.02 nm [7] . A reasonable 
estimate for the typical combined uncertainty of the 
relative spectral locations of the interference filter 
and the laser output wavelength would seem, there- 
fore, to be in the range of 0.01-0.04 nm, although 
larger errors in the filter CWL would increase this 
range. The range of 0.01-0.04 nm corresponds to ap- 
proximately 1—3 cm -1 when considering use of the 
tripled Nd:YAG laser. For the filters of ~0.25 nm 
width considered here, this uncertainty in the spectral 
location of the CWL of the filter translates into an 
uncertainty in the required temperature correction 
over a 100 K range of up to ~3% (refer to Fig. 4 in 
[16]). By comparison, as described above, the com- 
bined uncertainty of the relative location of the 
interference filter with respect to the laser output 
wavelength is determined by the uncertainty of the 
atmospheric layer temperature and the uncertainty 
of the fitting of the theoretical and measured curves 
as demonstrated in Figs. 3 and 4. Each of these uncer- 
tainties was estimated to not exceed 0. 1-0.2 cm -1 . 
This is approximately an order of magnitude lower 
uncertainty than that achieved by normal laboratory 
procedures. This lower uncertainty in the filter spec- 
tral location limits the uncertainty in the required 
temperature correction to ~0.25% at a given atmos- 
pheric pressure. Pressure changes of 25-50 hPa in 
the laboratory environment where the laser is housed 
do not significantly increase this uncertainty esti- 
mate. Also, it should be noted that this low uncer- 
tainty in the spectral location of the channel 
passband translates into high stability of calibration. 
With the ability to align the interference filter with 
the spectral passband with such accuracy, repeated 
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experiments as described here can help to maintain 
calibration even if temperature, mechanical disturb- 
ance or aging produces changes in the center of the 
passband [K)] . Repeated experiments of this type also 
help to randomize a component of the systematic 
error budget of the lidar which is advised for develop- 
ing time series of water vapor data for trend studies 
[28,29]. For a mobile system, such as the ALVICE 
laboratory, where measurements may occur at vary- 
ing altitudes or where mechanical disturbances may 
occur in transit, it is advisable to perform the tuning 
procedure described here at each location where 
measurements are taken. 

As stated previously, Fig. 2 indicates that, for the 
filter considered, selection of a filter center passband 
of 3649.5 cm -1 reduces the temperature sensitivity 
of F h to approximately 1% over the 100 K range 
considered. But this sacrifices between 4%— 12% of 
potential signal intensity when compared to tuning 
the interference filter to the peak of the F H curve. 
An application of Raman lidar focused on boundary 
layer meteorology or convective storm development, 
such as the IHOP experiment in 2002 [ 30 ] could 
benefit from placement of the interference filter at a 
center location of 3652 cm -1 . Whereas, an experiment 
focused on upper tropospheric ice supersaturation in 
the tropics or trend studies could well benefit from a 
placement of the center of the filter at 3653.2 cm -1 
since this placement of the filter would result in 
approximately 10% more water vapor signal strength 
for the cold upper tropospheric measurements as 
compared with those that would be taken if the filter 
were placed at 3649.5 cm -1 . With the techniques out- 
lined here, either of these peak spectral locations 
could be achieved with accuracy of ~0. 1 cm -1 ensur- 
ing that the temperature correction needed for either 
filter position could be accurately applied. 

Figure 5 demonstrates that knowledge of the 
atmospheric temperature profile to within ±5 K can 
limit the uncertainty of the applied temperature 
correction, due to uncertainties in atmospheric 
temperatures, to approximately 0.2%— 0.4%. When 
combined with the uncertainty in the temperature 
correction due to uncertainty in spectral filter loca- 
tion, the propagated uncertainty in the temperature 
correction becomes approximately 0.5%. If a system- 
atic bias of 0.5% can be considered insignificant for a 
given application, then an approximate temperature 
profile can be used for a sufficiently accurate appli- 
cation of the temperature dependence correction. If, 
for example, the temperature correction needs to be 
applied over a 24 h dataset, it may be necessary to 
use 2 or even 3 approximate temperature profiles 
during a 24 h period, where an interpolation can be 
done between the 2 or 3 profiles, to increase the 
accuracy of the correction. Such an approach can be 
used to limit the uncertainty of the temperature cor- 
rection to less than 0.5% even for data acquired over 
an extended period under day and night conditions. 

Correction of systematic bias is a required compo- 
nent of good metrological practice [ 31 ] requiring a 


sound estimate of the additional uncertainty in the 
final data product that is created by applying any 
corrections. Because it results in higher signal 
strength and allows for a more accurate spectral 
positioning of the interference filter, the case being 
made here is to tune the Raman lidar water vapor 
channel to a point that maximizes signal strength 
but also that creates a small temperature dependent 
bias. The bias can then be corrected, using the pro- 
cedures outlined here, where the additional uncer- 
tainty due to the correction is then included in the 
uncertainty budget of the data. Considering together 
an uncertainty in spectral location of the filter and 
an uncertainty in atmospheric temperature of 5 K as 
a function of altitude, the total uncertainty in the 
temperature correction determined here remained 
below 0.5%. This is the additional uncertainty due to 
the temperature correction. Given that Raman lidar 
random uncertainties are typically in the 5%— 10% 
range or greater and that calibration uncertainty 
is approximately 5% or greater [24,32], an increase 
of 0.5% in the total uncertainty budget would seem 
insignificant for most applications of Raman water 
vapor lidar data. If lower uncertainty in the temper- 
ature correction is needed, this can be achieved by a 
finer resolution fitting of lidar measurements to the 
theoretical F H curve (Figs. 3 and 4) and by using a 
more accurate profile of atmospheric temperature. 

7. Summary and Conclusions 

Previous efforts to calculate the temperature 
dependence of the Raman water vapor lidar spec- 
trum for use in Raman lidar profiling of water vapor 
have required very accurate characterization of both 
the output laser wavelength and spectral passband 
characteristics of the water vapor channel. Typical 
accuracies in the determination of these two quan- 
tities result in an uncertainty in the temperature 
dependence of the Raman water vapor lidar meas- 
urement of up to 3%. We have described here an 
experimental technique that evaluates the convolu- 
tion of the output laser wavelength and the spectral 
passband characteristics together with a resultant 
accuracy that greatly exceeds that of typical labora- 
tory results. The technique involves performing 
atmospheric profiling with a sequentially tilted 
interference filter to evaluate the curve of F H as a 
function of tilt angle. It is not necessary to know the 
effective index of the interference filter; it is only nec- 
essary to be able to discern a clear peak in the curve. 
The peak in this curve along with an approximate 
knowledge of atmospheric temperature permits the 
peak of the interference filter transmission curve to 
be determined with accuracy of ~0. 1 cm -1 . By posi- 
tioning the filter at the peak in this intensity curve, 
maximum transmitted intensity at that temperature 
is achieved although a correction is needed for the 
temperature dependence of the measurement. De- 
pending on whether the Raman lidar is optimized 
for lower tropospheric or upper tropospheric mea- 
surements, the magnitude of the temperature 
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correction ranges between 6%— 9% for a 100 K change 
in atmospheric temperatures. Under many circum- 
stances the range of temperatures from the surface 
to the tropopause will be less than 100 K indicating 
that the magnitude of the correction is less, as was 
the case illustrated in Fig. 5. It was also shown that 
only an approximate knowledge of the actual atmos- 
pheric temperature profile is needed to limit uncer- 
tainties in the correction to typically less than 0.5%. 
This increase in uncertainty is likely tolerable for 
any application of Raman water vapor lidar mea- 
surements. It is our hope that the procedure outlined 
here will ease the efforts to provide accurate narrow- 
band Raman water vapor lidar measurements. A 
Mathematica notebook that demonstrates the calcu- 
lations used here is available from the lead author. 

Appendix: A 

The filters used in this study possessed CWL approx- 
imately 0.05 nm longer than the peak of the water 
vapor Q-branch (e.g., see Fig. 3 in [4]). Because of 
this, it was possible to scan the full set of wavenum- 
bers shown in the figures in this paper. If, on the 
other hand, the CWL of the filter when positioned 
normal to the optical path is closer to the peak of the 
water vapor spectrum, it may not be possible to cre- 
ate a single-peaked F H curve as shown in Figs. 2—4. 
For example, consider the same Filterl from Fig. 1 
above but with the center of the filter, when normal 
to the optical path, at 3653 cm -1 . With such a filter, 
wavenumbers larger than 3653 cannot be reached 
through a tilt-tuning exercise since tilting an inter- 
ference filter from the normal position causes the 
transmitted wavelength to decrease. The type of 
F h curve that would result using such a filter is illus- 
trated in Fig. 6. The same fitting exercise described 
in the paper can be done to determine spectral 
location of the filter but the shapes of the curves will 
be quite different. Note that this effect does exist for 
experiments based on scanning the output of a 
spectrometer. 
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Fig. 6. F h versus wavenumber plot that results for a tilting 
experiment using a 0.26 nm filter with passband centered at 
3653 cm -1 when normal to the optical path. Note that no wave- 
numbers above 3653 can be reached through the filter tilting 
experiment. 
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